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Abstract. This paper reports results from a study of the reaction pp — } pK°E + at beam momenta of pteam 
= 2950, 3059, and 3200 MeV/c (excess energies of e = 126, 161, and 206 MeV). Total cross sections were 
OO _ determined for all energies; a set of differential cross sections (Dalitz plots; invariant mass spectra of all 

two-body subsystems; angular distributions of all final state particles; distributions in helicity and Jackson 
frames) are presented for e = 161 MeV. The total cross sections are proportional to the volume of available 
three-body phase-space indicating that the transition matrix element does not change significantly in this 
. range of excess energies. It is concluded from the differential data that the reaction proceeds dominantly 

via the 7V(1710)Pn and/or iV(1720)Pi 3 resonance(s); AT(1650)Sn and Z\(1600)P 3 3 could also contribute. 
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PACS. 13.75.Cs Nucleon-nucleon interactions - 13.75.-n Hadron-induced low- and intermediate-energy 
reactions and scattering (energy < 10 GeV) - 13.75.Ev Hyperon-nucleon interactions - 25.40.Ve Other 
reactions above meson production thresholds (energies > 400 MeV) 
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5~ ■ 1 Introduction effort has led to excitation functions measured at excess 

energies (e = \fs — (m p + ma + my)) from only a few 
The study of associated strangeness production in pro- MeV above the threshold to e « 250 MeV. Dalitz plots 
ton-proton collisions is one of the major physics programs were investigated by the COSY-TOF collaboration 
carried out at the COoler SYnchrotron COSY located at which, very recently, also published differential cross sec- 
Forschungszentrum Jiilich, Germany. Various experimen- tion data for both reaction channels [10] . Along with the 
tal groups have contributed data to the final states pK + A, wealth of data, various theoretical approaches based on 
pK+S°, and nK + £ + in the past decade [ UH I 3 l l l M 6 l |7 l very dlfferent footings were proposed [l^llTIl lTMTSl I 1711 181 
IM MTOlfni[l"2] . In all these three cases the final states con- [H]. In the case of the nK + S + final state total cross sec- 
tain a rather long-lived charged kaon {ct k + = 3.7 m) and tions were published for e < 200 MeV [7UHEI1II1] which, 
a nuclcon, and thus are experimentally well suited to be however, are contradicting each other strongly. This odd- 
accessed in inclusive and exclusive measurements at the ity is subject to an ongoing theoretical discussion [20] . 
COSY facility. In the case of A and S° production this 

Compared to the three final states discussed above the 

Correspondence to: m.schulte-wissermann@physik.tu-dresden.de experimental data basis for the reaction pp -» pK°S + is 

(M. Schulte-Wissermann) very scarce for e < 350 MeV. In fact, the only exper- 
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imental data for this reaction stem from the search for 
the supposed pentaquark state + (pp — > E + + ,0 + — » 
pK°) carried out by the COSY-TOF collaboration [21]. 
At higher excess energies measurements of total cross sec- 
tions were published originating from the bubble chamber 
experiments in the 1960s [^1^51[^l[^51[251[?7l[251E9] (com- 
piled in [30]). These data show, although with a large scat- 
ter, the total cross section to be rather constant (~ 20 fib) 
over a wide energy range (350 MeV < e < 2000 MeV). 
One of these early experiments reports a Dalitz plot [23] 
for e = 723 MeV. The approximately 30 entries show an 
enhanced density at small K°U + invariant masses which 
was interpreted as a hint for pion exchange. 

Also from the theoretical point of view, the reaction 
pp — > pK°U + was only poorly addressed as compared to 
the reaction channels containing a K + meson. In 1960, 
Ferrari [31] was the first to predict the total cross sec- 
tion for excess energies between 72 and 400 MeV in a me- 
son exchange model where pion as well as kaon exchange 
were considered. A schematic diagram of the respective 
exchange graphs is shown in fig. [T] on the left for kaon 
exchange (strangeness manifests itself in the exchanged 
boson) and in the middle for pion exchange (associated 
strangeness production at the pir — > K°S + vertex). An- 
other early calculation of the cross sections simultaneously 
for all pp — > NKY channels was performed by Tsu Yao 
[52] for a beam energy of T p — 2.85 GeV. This model 
is based solely on single pion exchange and the calcu- 
lated total cross section for pp — > pK°S + (e = 353 MeV) 
is in rather good agreement with the prediction of Fer- 
rari [31| and the only experimental value |22j known at 
that time. Later on, theoretical analyses and interpreta- 
tions appeared mainly in experimental papers [23 , 24, 26, 
127] . In 1968, Ferrari and Serio [55] explained fairly well 
in a meson exchange model all the then known total and 
differential cross sections for the NKY final state by in- 
troducing empirical cutoff factors in order to model var- 
ious damping effects at large momentum transfer (form 
factors, absorption). Since then, the theoretical progress 
in describing the production of the pK°S + final state 
was strongly hampered by the lack of new experimental 
data. In view of the experimental and theoretical progress 
made for the pK + A and pK + S° final states, in partic- 
ular data at lower excess energies are highly desirable. 
This would render possible the further development of 
more recent theoretical approaches 34,35,36 , which now 
also include, apart from kaon and pion exchange, a pro- 
duction scenario involving intermediate baryon resonances 
{pp -> B*p; B* = N*,A*; B* -> K°S+, see fig. [lb). 
If it is assumed that the associated strangeness produc- 
tion in pp-reactions procedes via resonances, a compre- 
hensive study needs the consideration of N* as well as 
A + * resonances in the case of both pp — > pK + U° and 
pp — > pK°U + reactions; due to isospin conservation the 
reaction pp pK + A can involve only N* resonances. 
The pp — » nK + Z" + -reaction, however, can proceed only 
via Z\ ++ *-resonances. 

In this paper we report results obtained for the reac- 
tion pp — > pK°U + measured for excess energies of e — 
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Fig. 1. Reaction mechanisms involving (a) strange, (b) non- 
resonant non-strange, and (c) resonant, non-strange meson ex- 
change for pp — >• pK°E + . Initial and final state interactions are 
not indicated. 



126, 161, 206 MeV (p beam = 2950,3059, and 3200MeV/c). 
These data supplement earlier studies on the reaction chan- 
nels pp —> ppuj [571158] . pp —> pK + A, and pp — > pK + S° 
|10j . As the data for these various reactions were taken 
simultaneously and the analyses utilized the same soft- 
ware package |3j5] the results obtained are characterized 
by a very high degree of internal consistency. This is, 
in particular, important for the future development of 
theoretic models when the three channels with associ- 
ated strangeness production are to be described simul- 
taneously. 

For all three excess energies total cross sections will be 
given. For the data taken at e = 161 MeV Dalitz plots, in- 
variant mass spectra, distributions in the CMS, in helicity 
as well as in Jackson frames will be presented. All differ- 
ential distributions are scrutinized with the aim of gaining 
insight into the reaction mechanism. It will be presented 
that all differential distributions can be described simul- 
taneously and consistently if a resonant production of the 
final state is assumed. 



2 Experimental procedure 
2.1 Detector setup 

The experiments were carried out with the time-of-fiight 
detector COSY-TOF located at an external beam line 
of the COoler SYnchrotron COSY (Forschungszentrum 
Jiilich). The COSY machine provides proton beams of 
very high quality (spill length 5 min; several 10 6 pro- 
tons/s; low emittance of < 5 7rmmmrad; relative momen- 
tum uncertainty Apjp < 10 -3 ). 

The layout of the COSY-TOF detector is shown in 
the upper part of fig. [2J in the lower part the near tar- 
get region with the time-of-fiight start and tracking detec- 
tors [4"UIl4"T] is sketched. The interaction volume is small 
and well defined as the narrow beam with Gaussian pro- 
file (cr x ,y < 300 /im) is directed onto a liquid hydrogen 
target of only 4 mm length [32] ■ The emerging particles 
traverse just behind the target (« 25 mm) a 24-fold seg- 
mented scintillation detector ( "start-detector" ) which pro- 
vides the start signal for the time-of-flight measurement. 
At a distance of 30 mm downstream of the target a double- 
sided silicon-microstrip detector is installed, followed by 
two double-layered scintillating fiber hodoscopes at 100 
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Fig. 2. The COSY-TOF detector (top), the near-target region 
(start-detector system, bottom). The lower picture shows a 
typical event pattern for the pK°E + exit channel. 



and 200 mm. These three tracking detectors measure the 
coordinates of traversing charged particles in three dimen- 
sions with a spatial resolution of ps 100 \im (microstrip) 
and w 1.5 mm (hodoscopes). 

After a flight path of ~ 3 m through the evacuated ves- 
sel (0.2 Pa) all charged particles are detected in the highly 
segmented stop components. They consist of two triple- 
layered forward hodoscopes (central and ring hodoscope) 
|%3"] and the barrel hodoscope [33] , all manufactured from 
BC412 scintillating material. From the combined measure- 
ment of time and position the velocity vectors of charged 
particles (originating from the target) are determined with 
a time-of-fhght resolution of better than cttof = 300 ps 
and an angular track-resolution of better than = 0.3°. 
Vertices from neutral particles decaying behind the mi- 
crostrip detector and before the first hodoscope are recon- 
structed from the tracks of their two charged daughter par- 



ticles with an accuracy of a XtV < 1 mm and a z < 3.0 mm. 
For the present analysis this feature is of crucial impor- 
tance as a measurement of the four-momentum of the Kg 
is mandatory for the full kinematic reconstruction of the 
final state. 

The COSY-TOF detector stands out for its low mass 
areal density of target, start-detector, and tracking detec- 
tors. This renders the influence of small angle scattering 
and energy loss almost negligible. In addition, the COSY- 
TOF detector has a high efficiency of > 95% for the de- 
tection of charged particles and covers a large solid angle 
(1° < < 60°, 0° < (j) < 360°) in the laboratory frame. 
These features allow the study of different reaction chan- 
nels (e.g. pp — > pp [?5], pp7 [3S], Wl [33, d7t + [33], ppw 
[35], pK°E+ |2B, and pK+A as well as pK+S° QD]) from 
the same data sample by examining the measured time-of- 
flight of the charged particles and the overall event topol- 
ogy- 



2.2 Principle of measurement and data analysis 

The TOF start-detector setup was designed to provide 
an effective means for the analysis of final states with 
open strangeness (e.g. pK + A, pK + U°, pK°E + ). Here, 
two charged particles are emitted at the three particle 
production vertex. In the case of pp —> pK°U + the K\ 
escapes the detector (ct k o = 15.3 m) while the K s has a 
considerable probability (ct k o = 26.8 mm) to decay into 
two charged particles behind the start-detector and be- 
fore the first hodoscope (see fig.[U lower part). Hence, the 
start-detector will be hit by only two, while the detectors 
located further downstream, i.e. the two hodoscopes and 
the stop detectors, will be hit by four charged particles. 
This "multiplicity jump of charged particles" is set as a 
trigger condition during data taking and is also the first 
condition required in the off-line analysis. 

In order to discriminate the reaction pp — > pK°S + 
from background, the characteristics of its final states is 
exploited: 1) a prompt track emerging from the target 
(proton), 2) a decay "V"(K® — > 7r + 7r~) with its vertex 
located in the decay volume (displaced vertex), and 3) 
one additional hit somewhere in the stop detector (due 
to the charged decay particle from either S + — > pir° or 
£+ -> nir+). 

Primary track and decay V candidates are selected by 
applying the following conditions: a primary track must 
have a signal in the start-detector and in one of the stop 
components. In addition, at least two signals from the 
three sub-detectors (microstrip detector, small and large 
hodoscope) are required. A secondary V must consist of 
two independent arms, each defined by fitting points in 
both hodoscopes and the stop detector. The point of clos- 
est approach of both arms is considered to be the kaon 
decay vertex; the vector connecting the center of the tar- 
get and this decay vertex is considered to be the kaon's 
direction of flight (j>k°)- I n order to discriminate a sec- 
ondary V candidate against the background generated by 
primary particles inducing reactions in the start-detector 
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a minimum angle of the pion with respect to the mother 
particle of 10° and a minimum angle between both pions 
of 30° is required. 

All permutations of primary tracks and decay V can- 
didates are subjected to an overall quality check which 
includes the number of involved fitting points, the dis- 
tance of closest approach of both tracks of the decay V, 
the quality of all fittings procedures (x 2 values), and the 
coplanarity of the decay V with respect to the primary 
kaon ($>k% x (f>7r+ x Pit-) ~ 0). The combination of a pri- 
mary track and a secondary V with the best overall quality 
is kept for further analyses. This method was developed 
by means of Monte Carlo data which show that in the 
final event sample 86% of the events are reconstructed 
correctly. 

So far only geometric information is exploited in order 
to identify the primary proton and the two decay pions 
of the primary (neutral) kaon. The velocity vector of a 
primary track is calculated from (t sto P — t start ^ anc i by as- 
signing the proton mass to this velocity vector the proton 
four-momentum is obtained. 

As the secondary pions from the K s decay emerge be- 
hind the start detector they lack an individual start-time 
information. Nevertheless, their momenta can be calcu- 
lated from the kinematic relations between the measured 
angels of the daughter pions with respect to their mother 
kaon. The sum of both pion four-momenta then yields the 
four-momentum of the kaon. The four momenta of beam, 
target, kaon, and proton are used in order to calculate the 
four momentum of the S + (Vs+ ^Vb + Vt — Vk° —'Pp)- 
The invariant mass of the S + is the "missing mass" of 
the final state. The resulting missing mass spectrum is 
the backbone of the analysis. 

Although the two pions do not provide an individual 
start-time signal their stop-time information can never- 
theless be used as a means to substantially reduce experi- 
mental background [3E] . For this purpose the difference of 
the time-of-flights of the two pions is calculated via: 

A tof = (C p - C rt ) - (C p - C rt ) (!) 

±stop istop 

as both pions are created simultaneously (t^ rt = t^ rt ). 
This measured time difference is then compared to the 
time difference determined from the path length of the 
pions in the detector and the pions' calculated momenta. 
Measured and calculated time differences must match with- 
in 2 ns {\A c t ff - A™f s \ < 2 ns). This requirement does 
not change the number of identified events significantly in 
the final event sample. The amount of experimental back- 
ground, however, is reduced by about a factor of two. 

It should be noted that the information of the E + 
enters only by the presence of a fourth charged parti- 
cle somewhere in the detector. The characteristic "kink- 
angle" (cf. lower part of fig. O could be used in order to 
discriminate between the two main decay channels of the 
hyperon S + — > pir° and E + — > nir + . In fact, due to the 
mass difference of the charged daughter particles 90% of 
the pir° (mr + ) decay branch is found with a kink-angle 



below (above) 12° . The ratio of the acceptance corrected 
counts is N s +^ p7T o /N s +^ nTr + = 0.99 ± 0.16, i.e. fully 
compatible with the accepted ratio of the branching ratios 
(B s +^ pn a /B s +^nn+ = 1-07 [35]). However, no selection 
based on the kink-angle has been applied in the final anal- 
ysis since it neither improves the missing-mass resolution 
nor increases the signal-to-background ratio substantially. 

The effects of various background reactions on the fi- 
nal missing mass spectrum were studied by Monte Carlo 
simulations. For this purpose 10 7 events, distributed ac- 
cording to equal population of the available phase-space, 
were generated for each of the reactions pp — > pp, dn + , 
ppTT , pnir + , ppuj, pK + S a and pK + A and subjected to 
the pp — ¥ pK°U + analysis routines. All but the hyperon 
channels produce negligible (< 1/10 6 ) or no background 
at all. The reactions pp — > pK + S° and pp — > pK + A cause 
a broad missing mass spectrum; no peak at the £ + mass is 
observed in either cases. The acceptance for both channels 
is about 1.5% with respect to the pK°2J + channel. Con- 
sidering the total cross sections for the A and E° channels 
[TO] their added contribution to the final event sample is 
below 15% for each of the three measured beam momenta. 

2.3 Acceptance correction and absolute normalization 

The Monte-Carlo package used [50,51 models the detec- 
tor and the physical processes to great detail. The event 
generator produces the particles of the exit channel either 
according to three-body phase space, or alternatively, in- 
termediate resonances can be chosen in order to model 
a two-step production process (pp — > pN* 7 N* — > K°E+). 
The particles (and their daughters, granddaughters, etc.) 
are then propagated through the detector. Branching ra- 
tios and lifetimes of all particles are incorporated accord- 
ing to the values given in [Jj5]. Energy loss, small-angle 
scattering, nuclear reactions, and 5-electrons are consid- 
ered. Digitized QCD- and TDC-signals are generated from 
the energy deposited in the active detector components. 
Noise and thresholds are modelled as known from the 
measured detector response. Deviations from an homoge- 
neously populated phase space can be introduced by weight 
functions on an event- by- event basis (this procedure is 
called filtering in the following) . Finally, the Monte Carlo 
data are subjected to the same routines as real data in 
order to determine the acceptanc^E 

The overall acceptance of about one percent for the 
reaction under study is determined by the ratio of posi- 
tively identified events to the number of generated (10 7 ) 
Monte Carlo events. This value is well explained by the 
event kinematics, the detector geometry, and the detec- 
tor performance: the probability to observe the neutral 
kaon via its Kg component (50%), the branching ratio 
of the kaon to two charged pions (69.2%), the probabil- 
ity of the secondary vertex to be located within the fidu- 
cial "decay volume" (31%), the probability of both sec- 
ondary pions to traverse the active detector range (23%), 

1 Throughout this paper, the term "acceptance" is used 
for the convolution of solid angle coverage, detector-, and 
reconstruction-efficiency. 
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Fig. 3. Proton- K° missing mass spectra measured for the three excess energies. A prominent peak is observed above a smooth 
background in all cases. The results of a fitting procedure for signal (dotted line), background (dashed line), and total spectrum 
(solid line) are shown in all cases. In the lower right frame the Monte Carlo result for e = 161 MeV is depicted. 



the reconstruction efficiency for both pions (60%), the re- 
construction efficiency of a primary proton track (95%), 
and the detection probability of the decay particle of the 
hyperon (82%). The acceptance has been subject to a de- 
tailed reevaluation, resulting in a modification of the total 
cross section at e=161 MeV presented below compared to 
the previously published value [21] . 

The relative uncertainty of the overall acceptance cor- 
rection was determined by considering the following ef- 
fects: 1) The efficiency for the detection of charged parti- 
cles is generally high (>95%) and known with an uncer- 
tainty of ±5% for each of the ss 1500 individual detector 
channels. This leads to a contribution to the uncertainty 
of the acceptance correction of 14% for the present re- 
action. 2) By altering all restrictions imposed during the 
data analysis the number of pK®£ + events in the final 
event sample can be changed by +10%/-70% while the 
total cross section changes by less than 10%. This value 
is taken as the contribution of the choice of the restric- 
tions used in the analysis to the overall uncertainty of the 
acceptance correction. 3) The influence of the choice of 
intermediate nucleon resonances (MC input) and angular 
distributions (filter) was found to be small (2% each), if 
the mass and width of the resonance and the asymmetry 
of the filter functions are altered within the limits imposed 



by the measurements (see below). Adding these contribu- 
tions quadratically an overall systematic uncertainty of 
18% is obtained. 



With regard to differential distributions the acceptance 
varies in all cases smoothly with the observable under con- 
sideration (details will be shown later when presenting the 
final results). Here, an additional uncertainty Adi comes 
about for each bin i due to the gradient of the acceptance 
function ai. This is accounted for by choosing as uncer- 
tainty the mean change of with regard to its adjacent 
bins, i.e. Z\a; = |(|aj - ai_i|/2 + \ai - a i+1 \/2). 



The absolute normalization is determined via the ana- 
lysis of elastic scattering, which was recorded simultane- 
ously during the experiment. Our results are normalized to 
the elastic scattering cross sections from the EDDA col- 
laboration |52j and yielded time-integrated luminosities 
of le.Qnb" 1 (e = 126 MeV), 214nb- 1 (161 MeV), and 
6.4 nb -1 (206 MeV). The total uncertainty of this proce- 
dure (5%) is due in equal parts to our analysis and the 
uncertainty of the literature data. For details see [3T?ll45| . 
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2.4 Determination of total and differential cross 
sections 

Figure |31 shows the pK° missing mass spectra obtained for 
the three excess energies of 127, 161, and 206 McV. Promi- 
nent signals for the E + hyperon can be seen at its cen- 
tral mass of 1189 MeV/c 2 above a smooth and structure- 
less background. As usual for time-of-flight detectors, the 
missing mass resolution (momentum resolution) is best 
for smaller velocities in the exit channel (smaller beam 
momenta in the entrance channel). When comparing the 
spectra one clearly notes the higher integrated luminos- 
ity for e = 161. At e = 126 and 206 MeV, however, a 
better beam quality was available, as manifested by the 
lower background contribution. In the lower right frame 
the Monte Carlo result is shown for e = 161 MeV. 

The number of events in the missing-mass peak are 
obtained by consecutively fitting first the background and 
then the signals. The background is parameterized by a 
quadratic polynomial, where only missing masses below 
and above the I7 + -peak are taken into account (m < 
1100MeV/c 2 and m > 1250MeV/c 2 ). Then, the back- 
ground parameters are fixed and the signal is described by 
a Voigt function (convolution of a Gauss- and a Lorentz- 
function). Voigt functions are chosen since they model 
properly the signal shape of a rather narrow peak accom- 
panied by broader tails. The integral of the Voigt function 
represents the total number of positively identified events. 
The overall systematic uncertainty due to signal and back- 
ground separation is determined by varying the fit-region 
for the background-fit below and above the S + -peak and 
was found to be ±4%. 

The total cross sections for all three excess energies 
are then obtained from the number of identified events 
(N), the integrated luminosity (£), and the acceptance 
(a) according to a — N/(C - a). As statistical uncertain- 
ties the numerical uncertainties of the fitting procedure of 
the signal will be quoted. As systematic uncertainties the 
quadratic sum of the uncertainties of luminosity determi- 
nation (5%), acceptance correction (18%), and signal-to- 
background separation (4%) will be given. 

Differential cross sections are determined in analogy 
to the total yield, only that the amount of signal and 
background is determined from individual missing mass 
spectra generated for each bin of the observable under 
study. Also the uncertainty of the signal-background sep- 
aration is determinded individually. Data will be shown 
if the uncertainty in a specific bin (root-mean-square of 
the statistical uncertainty combined with the differential 
uncertainty of the acceptance correction and the signal- 
background separation) is below 80% of its cross section 
value. 

The numerical values of all total and differential cross 
sections will be given below in tables. For the latter case 
the coefficients of least square fitting with Legendre poly- 
nomials 

da/dQ = ]T a r Pi, 1 = 0,1,2. (2) 



will be given in order to judge asymmetries (Pi, repre- 
sentative for all P dd) and anisotropies (P 2 ). It should be 
noted in passing that the integrals of all angular distri- 
butions (a tot = J j^df2 = 4n ■ do) match the total cross 
section within uncertainty. 



3 Results and discussion 

For a reaction of type a + b — » 1 + 2 + 3 the reaction space 
given by kinematics is 12 dimensional (three 4-momentum 
vectors of the exit channel), if the entrance channel is 
fixed (masses, y/s). Knowing the masses of the exit chan- 
nel reduces the dimensionality to 9 (three 3-momentum 
vectors). Four energy- momentum conservation equations 
reduce the dimensionality to five. If no spin direction is 
preferred, the azimuthal dependence is trivial and four di- 
mensions (variables) remain in order to uniquely describe 
the reaction kinematics. 

The actual choice of the linearly independent basis of 
such a four dimensional space (or four experimental ob- 
servables) is not unique. One convenient choice are inde- 
pendent Mandelstam-like invariants [53] : 

si = s 12 = (V 1 +V 2 ) 2 (3) 

s 2 = s 23 = (v 2 + r 3 ) 2 

t l =t a i = {V a -Vif 

t 2 = t b3 = (V b -V 3 ) 2 , 

where the Vi denote the four- momentum of particle i. The 
physical meaning of these variables can be infered from 
fig- HI ti, h are the squared momentum transfers between 
the initial particles and two of the three ejectiles; s\, s 2 
are the squared invariant masses of two of the three final 
state subsystems. In case of resonant production we adopt 
the convention that t\ is the momentum carried by the 
exchange particle while t 2 is the momentum transfer from 
particle b to one of the decay products of the resonance. In 
the case of a symmetric entrance channel as in the current 
case particles a and b are interchangeable. 




tj 



Fig. 4. Graphical representation of s, si, S2, ti, following the 
prescription of ref. [53] for the most general case of a ab — > 123 
reaction (left), and for the case of a two step reaction involving 
an intermediate resonance (right). 

Using these four variables the four-fold differential cross 
section can be written as 
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d 4 a/dsids 2 dtidt 2 = @ • \M(s; Si, s 2 , t\, h)f 



(4) 



where <P represents the properly normalized phase-space 
factor and M is the transition matrix element. 

The four relativistic invariants in cq. 0] are linearly 
connected to specific angles, for instance, s\ is a linear 
function of cos(Z(pi,p2)) in a properly chosen reference 
frame. Similarly, other angles can be used to substitute s 2 , 
t%, t 2 . Therefore, the four- fold differential cross section can 
also be written as a function of a combinations of angles 
and invariants. 

The experimental goal is to provide the four-fold dif- 
ferential cross section (in whatever basis). This, however, 
is unrealistic in many cases, simply for statistical reasons. 
In practice, one therefore needs to reduce the dimension- 
ality by projecting onto subspaces: Integration over t\ and 
t 2 results in Dalitz plots. Projecting onto one dimension 
yields, for instance, invariant masses or angular distribu- 
tions in particular reference frames. 

The reduction of dimensionality is accompanied with a 
significant loss of information, as possible correlations can 
no longer be recognized. In addition, one has to be cau- 
tious about the kinematic correlation of the three body 
final state, as in specific cases a true physical cause in 
one variable can mimic a characteristic signal in another 
variable [53] (kinematic reflection). It therefore is essential 
to evaluate as many as possible (linear independent) pro- 
jections of the four-fold differential cross section. These 
projections are then to be described simultaneously by a 
theoretic model. 



3.1 Total cross sections 

Integration over all four variables in eq. (j4j) results in the 
total cross section. Thus, a single value for the total cross 
section does not allow to infer any details of the reac- 
tion mechanism. Nevertheless, the evolution of the total 
cross section with excess energy, the excitation function 
(<7 = a{yfs)), is often used as a first means to judge how 
well (different) theoretical approaches are appropriate to 
describe these data. 

The results for the total cross section are listed in table 
[TJ They are included in fig. [5] which show the world data 
for the reaction pp — > pK°E + . All data above e = 300 
MeV stem from bubble chamber experiments [22,23,24, 
I25U26U27U28LI29U30J . The data show a relatively large scat- 
ter, however indicate a rather constant total cross section 
over a wide energy range. The present data are the first 
to establish the excitation function at smaller excess en- 
ergies. 

The results of two early theoretical approaches are 
shown in the same figure. In 1962 Tsu Yao [33] calcu- 
lated the total cross section for the present reaction at 
e = 353 MeV considering solely single pion exchange. The 
resulting value of a — 36 pb is in good agreement with the 
only experimental value [35] available at that time. The 
calculation by Ferrari and Serio in 1968 [33 (based on the 



Table 1. Total cross sections for the reactions pp — ► pK°E + 
for the three different excess energies. The first uncertainty 
refers to statistical and the second to systematical ones. 



e (MeV) 


acceptance (%) 


signal (#) 


o-tot (/ib) 


127 


0.93 


386 


2.46 ± 0.13 ± 0.47 


161 


1.05 


9226 


4.13 ± 0.06 ± 0.79 


206 


0.93 


412 


7.02 ± 0.36 ± 1.34 
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this work 

bubble chamber data 
Tsu Yao 1962 
Ferrari & Serio 1968 
Tsushima et al. 1999 
phase space 
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Fig. 5. Total cross sections of the reaction pp — » pK°E + . The 
triangles at low excess energies represent the present data. The 
solid square symbols are from bubble chamber experiments 
(compiled in 30 ). The theoretical calculations shown are de- 
scribed in the text. 



pioneering work of Ferrari |31j ) was aimed at the simulta- 
neous description of various pp — > NKY channels within 
a boson exchange model (explicitly excluding nucleon res- 
onances) . The energy dependence of the total cross section 
calculated by this model is shown as dotted lines for excess 
energies up to 1500 MeV. Here, the lower line represents 
the results for purely pion exchange, while the additional 
contribution of kaon exchange (for a coupling constant 
of G\ less than or equal to 1.6) is represented by the 
hatched area. The calculation for pure pion exchange is in 
rather good agreement with our data. Nevertheless, it will 
be shown in the following that nucleon resonances play 
an important role for the production of the pK°S + final 
state. Hence, the model of Ferrari and Serio is lacking a 
central ingredient and the good agreement between the 
experimental data and theoretical prediction is probably 
a coincidence. 

More recent calculations by Tsushima et al. [35] are 
based on a resonance model where it is assumed that 
an excited baryon, B* (N* or A*), is produced via me- 
son (tt, r),p) exchange (pp — > pB*; B* — > K + A, K + U a , 
K°S + ). With one set of parameters (coupling constants, 
cut-off values) all these three hyperon channels are treated 
on equal footing. In the case of pp — > pK°S + this model 
describes the total cross sections measured at excess ener- 
gies above 800 MeV rather well (dashed curve), however 
misses the experimental value at e =353 MeV by about a 
factor four. This model underestimates the new data by a 
factor two. 
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m| v (GeV 2 /c 4 ) m 2 j,. (GeV 2 /c 4 ) 



Fig. 6. Dalitz plots of the squared invariant masses of the 
K°p vs K°E + (left) and K°E + vs pE+ (right) subsystems for 
e = 161 MeV. The data is acceptance corrected, however the 
background is not subtracted (see text). In both frames a clear 
enhancement at low K° S + masses indicate the presence of an 
intermediate resonance state (B* — > K°E + ) located near the 
lower mass boundary. A contribution of a pE + FSI would lead 
to an enhancement of data at the lower mass boundary of the 
pS + mass axis. A significant signal of pS + FSI is not observed. 



sons. However, sideband cuts show that the structures do 
not originate from the background. 

In both Dalitz plots structures due to narrow reso- 
nances or final-state-interaction are not observed within 
resolution. However, the relative bin occupancy of both 
Dalitz plots continuously increases towards low K°E + - 
masses. This can be interpreted to be caused by one or 
more intermediate resonances (B* — » K°£ + ) with cen- 
tral masses (m;) in the region of the lower mass boundary 
(or below) and broad widths (P,) of at least 100 MeV. The 
PDG [35] offers a whole list of baryon resonances which 
could potentially interfere and contribute to the reaction. 
Due to the arguments given above, it is not trivial to ex- 
tract the resonance parameters (to^-Z^) from the Dalitz 
plots alone. In the following, however, we are going to ar- 
gue that for the specific case of pp — > pK°S + a set of 
one dimensional distributions is well suited to shed light 
on the reaction mechanism, particularly on the resonant 
contribution. 



We also include a solid line which represents the vol- 
ume of available phase space (a oc a ■ e 2 , a = 1.617 • 
10~ 4 nb MeV -1 ). It describes the data very well within 
the given experimental uncertainty. This indicates that 
the absolute value of the transition matrix element in eq. U 
does not depend strongly on excess energy in the inves- 
tigated energy region. However, this does not imply that 
the transition matrix element itself is constant over the 
available phase space. 



3.2 Dalitz plot 

Integrating over t\ and ti in eq. [4] results in the double- 
differential cross section cPa / 'ds\ds2, whose representation 
is the Dalitz plot It connects two Lorentz invariants 
{e.g. si = m\ 2 , s 2 = m 23i squared invariant masses) 
and, hence, gives insight into the correlation of the three 
particle final state. If the matrix element in eq. 0] is con- 
stant, i.e. does not depend on s\, S2, t\, t 2 , the Dalitz 
plot is homogeneously populated. If, however, the matrix 
element is not constant, the deviation from an uniformly 
populated phase-space may show up as characteristic dis- 
tortions. Prominent examples are resonances and FSI ef- 
fects. If these signals are well separated, direct conclusions 
on physical properties can be drawn (e.g. mass and width 
of a resonance; scattering lengths via FSI). This simple in- 
terpretation is not possible if different effects overlap on 
the Dalitz plot. In this case a proper theoretical model 
(treating the dependence of the matrix element on all four 
variables s±, s 2 , ii , t 2 ) must consider the different contri- 
butions coherently on the amplitude level. 

Acceptance corrected Dalitz plots (w^Op vs m^ 0s+ 

and Wko£+ vs TO pS+) are snown m fig- [6] for an excess en- 
ergy of 161 MeV. Data is shown for a missing mass region 
of ±40MeV/c 2 around the S + mass, thereby reducing the 
contribution of background to below 20%. A bin- wise sub- 
straction of background is not feasible for statistical rea- 



3.3 Angular distributions in the helicity-frames and 
invariant mass distributions 

Invariant-mass distributions are obtained by projecting 
the Dalitz plot onto one of its axes, s\ = m 2 2 , S2 = m 23i 
or S3 = m 2 ^. Although the resonance parameters (m,r) 
are reflected in the invariant mass spectra, their direct 
extraction is only possible for an isolated non interfering 
resonance. In any case, one has to be cautious about kine- 
matic reflections as the three invariants are related via 
si + s 2 + S3 = s + to 2 + m 2 + m\. Any structure deviating 
from phase space in e.g. s 2 , will cause deviations of the 
pure phase-space behavior of the other observables si and 
S3. It is therefore not a priori clear which invariant shows 
a cause and which an effect. In contrast, the situation is 
more clear when analyzing the angular distributions in 
helicity frames. 

Angular distributions in helicity framed are essentially 
special projections of a Dalitz plot. There is a linear re- 
lation between Si and cos9^ k (and cyclic permutations), 
therefore si and S2 can be substituted by two helicity an- 
gles in eq. 0] [53] . 

For a single resonance (in the following we assume it 
to decay into the 23-system) the properties of mass and 
width lead to anisotropic distributions in the 12- and 13- 
hclicity frames, while the angular distribution in the 23- 
helicity frame remains isotropic (projection of the "reso- 
nance band" in the Dalitz plot). This isotropy in the 23- 
hclicity frame is independent of mass and width of the res- 
onance and therefore uniquely identifies the decay chan- 
nel. This statement is also valid if more than one resonance 



2 For reactions of type ab — s- 123 the 23-helicity-frame (R23, 
R indicates .Reference frame) is defined as the rest frame of 
the particles (2,3) (ps = — P2). The respective (polar) helicity 
angle in this frame is defined as the angle between particle 3 
and 1 (8\i 3 ). By cyclic permutation three helicity frames can 
be constructed for the three-body final state (R23, R31, R12). 
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Fig. 7. Angular distributions of the particles in helicity frames (top) and the invariant mass distributions for all three two- 
body subsystems (bottom) measured at an excess energy of e = 161. Error bars for each data point are the square root of the 
quadratic sum of the statistical, acceptance, and signal-to-background-separation uncertainty. The dashed histograms represent 
the effect of an N(1720,150) intermediate resonance, which is used as input for the Monte Carlo simulation. The grey curves 
show the distibutions expected for a pure phase-space scenario. Below each distribution the differential acceptance is shown. 
The numerical values of the differential cross sections are listed in table U and 



(decaying into the same two-body system) contribute in- 
coherently. 

The angular distributions in all three helicity frames 
measured at an excess energy of e = 161 MeV are shown in 
the upper row of fig. [71 In the lower row the spectra of all 
three invariant-mass subsystems are displayed. The accep- 
tances shown under each distribution are rather constant 
for all spectra. Due to the the method of determining the 
yield the physical background is subtracted individually 
for each data bin (cf. sec. 12. 4[) . 

It is obvious that the distributions in the pS + and 
K°p helicity frames are asymmetric (i.e. a\ ^ 0) while 
that in the S + K° frame is almost isotropic (ai ~ 0). 
This can be infered quantitatively from the values of the 
Legendre polynomials coefficients listed in table [5J It is 
also evident that the K°U + invariant mass distribution 
deviates from that given by phase space (grey curve in 
fig. 7) and shows a clear enhancement towards smaller 
masses; thus, strongly advocating a production scenario 
involving a broad intermediate resonance with a central 
mass near (or below) the lower K°S + mass boundary. 

The isotropic distribution observed in the K°S + he- 
licity frame is of particular relevance in case of resonant 
production: Firstly it indicates that the decay channel of 
the resonance is B* — > K°S + . Secondly, there is no pre- 
ferred orientation of the £?*-spin with respect to the B* 
direction of flight. This is an indication of the presence of 
several partial waves. And thirdly, the observed isotropy 
signifies that only one resonance participates or, if more 
than one resonance is involved, they do not interfere (ex- 



cept for the possibility that the interference pattern mim- 
ics an isotropic distribution). 

In the following, we assume that only one resonance 
participates. Indeed, all six distributions of fig. [7] are well 
described by the results of a Monte Carlo simulation (dashed 
histograms), where the three-particle final-state is mod- 
eled via pp — ► pN*, N* -> K°E + . The central mass and the 
Breit-Wigner width of the resonance were chosen to be 
m N , = 1720MeV/c 2 and r = 150MeV/c 2 (abbreviated 
in the following as N(1720,150)). 

The choice of this particular N(1720,150) resulted from 
a series of Monte Carlo simulations which were performed 
to study the kinematical effect of mass and width of dif- 
ferent resonances such as N(1400,270) Q£], N(1650,300), 
N(1720,150), and N(1900,300) chosen from gH]. We also 
included an N(1800,200) state in order to study the devel- 
opment of the kinematical effect with a narrower spacing 
of masses. (As kinemantics does not depend on isospin, 
the following statements apply also for Z\(1600), 4(1620), 
4(1700), 4(1750), and 4(1900) g5].) 

As an pictorial example of such an analysis the distri- 
butions obtained in the pS + helicity helicity are shown in 
fig. [5] Due to kinematics, all heavy resonances have a pos- 
itive slope in the pS + helicity frame. They clearly fail to 
describe the data. In contrast, all distributions resulting 
from the decay of light resonances display negative slopes 
in this helicity frame and describe the data with similar 
accuracy. 

In order to assess the quality of the overall description 
of the data by each of these N* resonances Monte Carlo 
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Table 2. Legendre polynomial coefficients (in units of nb/sr) 
determined by least square fitting to angular distributions of 
the reaction pp -»• pK°S + at e = 161 MeV, in the overall CMS, 
the Jackson and helicity frames (top to bottom). 



cos 


ao 






a 2 


n* 
u p 


335 ± 20 


O ± 1 


oy ± oy 


9k 


330 ± 22 


-10 ±44 


109 ± 39 


U S+ 


325 ± 11 


22 ±20 


46 ±29 


nRKp 
VbK 


319 ± 


6 


38 ± 12 


52 ± 16 


nRpS+ 
% 


322 ± 


7 


-27 ±14 


76 ±20 


nRS+K 
d bE+ 


322 ± 


7 


-7 ±13 


99 ± 16 


aRKp 


323 ± 


9 


183 ± 16 


-20 ± 20 


oRpS+ 


326 ± 


9 


-130 ±17 


41 ±22 


aRS + K 
°pE+ 


332 ± 


9 


19 ±17 


-19 ±22 



simulations a simple quality criterion - root-mean-square 
of the error-weighted difference of data and Monte Carlo 
results, averaged over all bins for the six observables of 
fig. [7J- was applied. It turns out that this figure of merit 
varies by about 10% for the three light resonances, where 
the N(1720,150) yields the best averaged description. This 
resonance is used throughout the whole analysis as Monte 
Carlo input for acceptance correction. 

The working hypothesis of a single resonance governing 
the pK°U + reaction process is further applied in the next 
subsection. 
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Fig. 8. Angular distribution measured in the pU + helicity 
frame. The lines are the Monte Carlo result for various res- 
onances (N(m, r)). Simulations with low mass resonances re- 
produce the measured data with similar quality, whereas heavy 
resonances clearly fail. The resonance with the best overall 
agreement between data and Monte Carlo is the N(1720,150) 
which is plotted as solid line. 



3.4 Angular distributions in the CM- and 
Jackson-frames 

Angular distributions in the overall CMS can be described 
by Legendre polynomials. In the specific picture of a pro- 
duction mechanism involving an intermediate resonance 
(pp — > pB*) the CMS scattering angle of the proton di- 
rectly reflects the angular momenta involved in the pB* 
channel. There is a linear dependence of the cosine of the 
scattering angle on momentum transfer t\, thus the pro- 
ton CMS angle can be taken as one of the linear indepen- 
dent variables in eq. 2] Due to the symmetric pp entrance 
channel all CMS distributions must of course be symme- 
tric with respect to cos6>* = 0. 

Jackson framed are the natural Lorentz invariant frames 
to investigate the angular momenta involved in a two- 
particle subsystem [56 . To illustrate this, consider the 
secondary decay (B* -> K°S + ) at the 7rp -> K°U + ver- 
tex in fig. [TJd, c; this represents a "2 — > 2" reaction. In 
this picture, the inverse reaction (K° S + — > pir) must have 
the same properties due to time reversal invariance. If one 
now imagines colliding beams of kaons and hyperons with 
Py = — Pk, which is by definition the KY Jackson frame, 
it is self-evident that the distribution of angles 0^ Y of 
the (in this case emerging) proton with respect to the (in 
this case beam-axis defining) kaon contains information on 
the relative angular momenta involved. There is a linear 
relation between and cos0?^ + s , i.e. the K°S + Jack- 
son angle can substitute the momentum transfer between 
the beam particle and one of the resonance decay parti- 
cles in eq. |U It is important to notice that the Jackson 
frame is a different Lorentz frame than the CMS and no 
symmetry with respect to cos# = is required. However, 
the distributions measured with respect to the target and 
beam proton must be identical. 

The angular distributions of all primary particles in 
the CMS measured at an excess energy of e = 161 MeV 
are shown in the upper row of fig. |H1 The lower row shows 
the angular distributions in all three Jackson frames. The 
variation of the acceptance is shown under each distribu- 
tion. In case of the proton and kaon the acceptance shows 
a quite strong angular dependence which is due to the 
strongly decreasing probability of the kaon to reach the 
decay volume when emitted further into the backward 
CMS hemisphere (cos8* K0 —> —1). This experimental ef- 
fect is mirrored in the proton acceptance distribution, as 
protons and kaon tend to be emitted into opposite hemi- 
spheres. For the hyperon, constructed from the combined 
momentum vectors of kaon and proton, a much weaker an- 



3 For reactions of type ab — > 123 the Jackson frame is de- 
fined as the Lorentz frame in which the center of mass of the 
particles (2,3) is at rest (p3 = — pz). The 23-Jackson frame and 
the 23-helicity frame are in fact the same Lorentz frame (_R23). 
The word "frame" rather refers to the choice of the reference 
axis which in case of Jackson frames is the direction of the 
beam particle. The Jackson angle is then defined as the angle 
between the beam direction and that of particle 3 (0^ 3 ). By 
cyclic permutation three Jackson frames can be constructed 
for the three-body final state (RES, R31, R12). 
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Fig. 9. Angular distributions of the particles in the overall CMS (top) and the Jackson frames (bottom), for the reaction 
pp — ► pK°E + measured at an excess energy of e = 161 MeV (pbeam = 3059MeV/c). Error bars for each data point are the 
square root of the quadratic sum of the statistical, acceptance, and signal-to-background-separation uncertainty. The solid 
histogram in the proton CMS and E + K° Jackson frame represents the respective Legendre polynomials of table [2] which are 
used as MC filter. Their affect on all other angular distributions are shown by the dashed histograms. Below each angular 
distributions the differential acceptance is shown. The numerical values of the differential cross sections are listed in table [S] 



gular dependence of its acceptance is found. Nevertheless, 
after acceptance correction all spectra show the necessary 
symmetry with respect to cos6** = (<zi = 0), as the val- 
ues of <zi given in table [5] are all compatible with zero 
within uncertainty. The e^/arj ratio changes by less than 
7% if ai is forced to be zero. 

A clear anisotropy is observed for all Jackson-frame 
distributions, as indicated by the a\ coefficients listed in 
table [5] There is no need to introduce Legendre polyno- 
mials of order higher than two. In the picture of a two 
step production process, this finding for the proton CMS 
angular distribution and the distribution in the K a S + 
Jackson frame indicates that angular momenta of I < 1 
are involved at both the pp -» pB* and B* — > K°E + 
subprocess. 

In order to corroborate the two step production sce- 
nario we now turn to the simultaneous description of the 
data by Monte Carlo simulations, the result of which is 
shown as either solid or dashed histograms in fig. [5] The 
solid histograms shown represent Monte Carlo data which 
were filtered by means of weight functions in the two natu- 
ral frames described above: The first filter was tailored to 
reproduce the measured proton CMS distribution (solid 
line in fig. [SI upper left frame). This filter significantly 
affects via kinematic correlation the two Jackson frames 
containing the proton (distribution in the pK° and pS + 
Jackson frames, dashed line in fig. [5J lower left and lower 
middle frame), while barely affecting the kaon and hy- 
peron CMS distributions. The angular distribution in the 
K°£ + Jackson frame is not affected at all. Hence it is 



justified to apply a second filter which is tailored to re- 
produce the measured distribution in the S + K° Jackson 
frame (solid line in fig. HO lower right frame). This filter sig- 
nificantly affects via kinematic correlations the kaon and 
hyperon CMS distributions (dashed lines in fig. [9j upper- 
middle and upper-right frame) while barely affecting the 
distribution in the pS + and pK° Jackson frames. The pro- 
ton angular distribution in the CMS is not affected at all. 
The observation that both filters do not effect each other 
is a consequence of the linear independence of t\ and . 

Two other filter combinations also fulfill the require- 
ment of linear independence (i^°-CMS and pS + Jackson 
frame; Z' + -CMS and pK°- Jackson frame). Although these 
filter combinations represent somewhat awkward reaction 
scenarios (i.e. pp —> K°X, X —> pS + and pp — s- S + X, 
X — > pK°) they are not excluded a priori. Therefore both 
combinations were tested for the sake of completeness with 
the method introduced in sec. 13.31 Both filter combina- 
tions yield a description of the data slightly inferior to 
the filter combination discussed above (the ones shown as 
solid histograms in fig. [9]) . 

In the following it is assumed that both the distribu- 
tion in the K°£ + Jackson frame and that of the proton in 
the overall CMS reflect the physical cause and the other 
four distributions are simply their kinematical reflections. 
In addition the discussion will be based on the proper- 
ties of resonant K °E + production from an initial state of 
identical particles. The angular distribution in the Jackson 
frame reveals that the K°S + subsystem carries angular 
momentum of I < 1 (cf. tabled]). 
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Table 3. Possible transitions for pp — ¥ pN(Sn), pN(Pu), pN(Pis) which yield differential cross sections described by da/dfl — 
fflo + a,2P2{cos9) with all a 2 coefficients being 7^ but those of transition la and 3e which have 02 = 0. Note that initial partial 
waves of type 1 So and 3 Po result in isotropic angular distributions irrespectively of the intermediate resonance. 
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Fig. 10. Selected theoretical angular distributions calculated 
for the pp — > pN* reaction (TV* : Sn, Pn, P13) according to 
|57] . All distributions are of the type da/dil = ao + a2p2(cosS), 
i.e. l P N* < 1. 



It then follows from isospin, angular momentum, and 
parity conservation that only nucleon resonances of the 
type S\i, P\i, and Pi 3 and Z\-resonances of the type S31, 
P31, and P33 can be involved. All other resonances have 
decay angular momenta of I > 1. 

In the next step, all angular distributions for protons in 
the CMS are calculated for pp pN(Sn), pN(Pn), and 
pN(Pi3) applying the formalism of Blatt and Biedenharn 
[57] . Entrance channel partial waves up to 3 Pj are consid- 
ered (spectroscopic notation 2S+1 Lj, with S: channel spin, 
L: orbital angular momentum of the two-body system, J: 
total angular momentum of the channel). Angular distri- 
butions arising from 1 So and 3 Pq entrance channel partial 
waves are isotropic irrespective of the final state due to 
J = for the entrance channel. 

All other transitions are listed in table [3J Interestingly, 
all angular distributions resulting from these transitions 
are of the type da/dft = ciq + a2P2(cosdcM), i-e. no term 
higher than Pi is involved which is perfectly in line with 
the experimental finding of I < 1. The calculations ac- 
cording to ref. [57] yield ai = for transition la and 3e 
(isotropy), a 2 > (concave shape) for lc, 2b, 2c, 2d, 3a, 
3c, 3f and a-i < (convex shape) for lb, 2a, 3b, 3d, 3g. 

In the case of Z\-resonances isospin selectivity requires 
their production from a spin singlet state (isospin triplet 



state) in the entrance channel, thus ruling out entrance 
channel P- or P-wave contributions. It is found that only 
two transitions are possible which yield angular distribu- 
tions with / < 1 namely { 1 D2)i — > ( 3 p2)/ involving the 
Z\(S3i) intermediate state and { 1 D 2 )i — > ( 5 S I 2)/ via the 
^(P33)j the latter being isotropic. A Z\(P 3 i) resonance 
yields angular distributions with an c^Pt term indicative 
of I — 2 which is not observed; hence such a resonance is 
excluded from the further discussion. 

Examples of calculated distributions resulting from tran- 
sitions listed in table [3] are shown in fig. [TU] normalized 
to the experimental data. The concave curves lc and 2b 
match the data best. However, it is very unlikely that only 
one or two specific transitions govern the reaction process. 
In fact, an inspection of table [3] reveals that the same ini- 
tial and final state partial waves can occur, however with 
different intermediate resonances being involved (e.g. 2c - 
3c, 2d - 3f). In addition, various initial state partial waves 
pass through a particular resonance but end up in differ- 
ent final state partial waves (e.g. 2a-2b, 3c-3d). A rather 
complicated scenario turns up and it is therefore impossi- 
ble to pin down the individual contributions of particular 
transitions within this simple reaction model. The ques- 
tion of which partial waves dominate the reaction can only 
be answered by a complete partial wave analysis. 

Nevertheless, fig. HU1 indicates that a proper superposi- 
tion of transitions will match the experiment data. Thus, 
the working hypothesis introduced in Sect. 13.31 of a re- 
action mechanism involving an intermediate resonance is 
now limited to the states Sn, Pn, P13, S31, and P33. 

The clear non-zero angular momentum observed in the 
K°S + Jackson frame requires considerable strength of 
participating resonances with intrinsic angular momen- 
tum, i.e. favoring Pn and/or P13 and/or P33. Resonances 
with these spins and parities can be related to those listed 
by the PDG. It turns out that 7V(1710)Pn, N(1720)Pi 3 , 
and Z\(1600)P33 are the only candidates with appropri- 
ate mass and a KS decay branch. Exactly this type of 
resonance (m ss 1720 MeV, P w 150 MeV or m « 1650 
MeV, P « 300 MeV) was also favored in describing the 
invariant mass and helicity angular distributions. As no 
Z\(S*3i)-resonance is known in this mass range which de- 
cays into K£, the ^(lGSO)^!! is the only further inter- 
mediate resonance to be possibly involved. 

It must be emphazised that the results for mass and 
width of a potential intermediate resonance obtained in 
sect. 13.31 and in this section resulted from exploiting or- 
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thogonal parameter spaces, namely s±, 82, s 3 and t\, 
respectively. The possibility of cross-checking results ob- 
tained in various reference frames and their combined in- 
terpretation underlines the importance of an analysis of 
the full four-dimensional reaction space. 



4 Summary 

The reaction pp — > pK°2J + was studied at the excess 
energies of e = 126, 161 and 206 MeV with the COSY- 
TOF spectrometer. The large acceptance allows the four- 
fold differential cross section d 4 a / ' ds\ds2dt\dt-2 to be de- 
termined. This high dimensionality however, needs to be 
reduced by projecting onto a reduced number of observ- 
ables to highlight certain aspects of this four-dimensional 
parameter space. In doing so, Dalitz plots, invariant mass 
spectra, angular distributions in helicity frames, the over- 
all CMS, and Jackson frames are obtained for e = 161 
MeV. 

The measured total cross sections are the first to build 
the excitation function between the production threshold 
and the energy region covered by early bubble chamber ex- 
periments. The observed increase of the total cross section 
with excess energy is consistent with the increase of phase 
space volume. This, however, does not mean that the re- 
action is governed by phase space; it rather implies that 
the transition matrix element does not (strongly) depend 
on energy up to about ss 200 MeV. The bubble chamber 
data show the total cross section to be rather constant 
(~ 20 pb) in the energy range from 300 to 1500 MeV. 

The Dalitz plots do not show any sign of final-state- 
interactions. However, a clear signal for resonant contri- 
butions {pp pB*; B* = N*,A*; B* -> K°U+) is 
found near the lower K°S + mass boundary. This qualita- 
tive result is confirmed by the quantitative analysis of in- 
variant mass spectra and angular distributions in helicity 
frames. The best description of the data by Monte Carlo 
is found for resonances with mjv* ~ 1720MeV/c 2 and 
P s» 150McV/c 2 . However, a contribution of lighter N* or 
A* resonances cannot be excluded. In contrast, the data 
clearly rule out a significant contribution of resonances 
with masses larger than 1800 MeV/c 2 . The isotropic dis- 
tribution in the K°E + helicity frame shows again the ab- 
sence of pS + final-state interactions as well as the ab- 
sence of interference effects which would distort this dis- 
tribution. The absence of interference allows the produc- 
tion mechanism to be interpreted as a two step process 
pp — > pB* , B* — > K° S + involving only one resonance. 

The decay angular momentum of the resonance is re- 
flected in the K°S + Jackson frame angular distribution 
and was found to be I < 1. Parity and angular momen- 
tum conservation then restrict the possible states to be 
Su, Pu, P13 (N*) or P 33 (A*). The shape of proton an- 
gular distributions in the CM system for the pp — > pB* 
reaction were calculated under the assumption that one of 
these resonances is involved. It turned out that all these 
calculated distributions have a relative angular momen- 
tum in the p — B* system of I < 1 in agreement with 



that deduced from the experimental proton CMS distri- 
bution. Considering that the angular distribution in the 
K°S + Jackson frame is clearly anisotropic and in view of 
the conclusions drawn from inspecting the distributions of 
invariant masses and helicity angles, JV(1710)Pn and/or 
A^(1720)Pi3 are the remaining candidate(s) for the dom- 
inating process of the pp — > pK°S + reaction. Contribu- 
tions from AT(1650)5u and/or Zi(1600)P 33 could also be 
present. 

A simultaneous description of all 12 differential dis- 
tributions measured was possible by assuming a single 
resonance (iV(1720), P = 150MeV/c 2 ) to dominate the 
reaction and applying weight functions (filters) on both 
the proton CMS angular distribution and the distribution 
in the K°S + Jackson frame. These kinematic constraints 
are linearly independent in the four-dimensional reaction 
space. Hence, this simple approach to describe the reac- 
tion mechanism is self-consistent. 

An advanced analysis of the data could be based on a 
partial wave analysis as performed by the Bonn-Gatchina 
group for various other reactions [58]. This might stimu- 
late the further development of theoretical models which 
describe the associated strangeness production in proton- 
proton collisions from first principles. 
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Table 4. Differential cross sections (in units of nb/(9.0MeV/c 2 )) of the invariant mass distributions shown in fig. [7] 
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Table 5. Differential cross sections (in units of nb/sr) of the angular distributions shown in fig. [7] and 
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